Reentrant phase diagram and pH effects in cross-linked gelatin gels 
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Experimental results have shown that the kinetics of bond formation in chemical crosslinking 
of gelatin solutions are strongly affected not only by gelatin and reactant concentrations but also 
by the solution pH. We present an extended numerical investigation of the phase diagram and of 
the kinetics of bond formation as a function of the pH, via Monte Carlo simulations of a lattice 
model for gelatin chains and reactant agent in solution. We find a reentrant phase diagram, namely 
gelation can be hindered either by loop formation at low reactant concentrations, or by saturation 
of chain active sites via formation of single bonds with a crosslinker at high reactant concentrations. 
The ratio of the characteristic times for the formation of the first and the second bond between 
the crosslinker and an active site of a chain is found to depend on the reactant reactivity, in good 
agreement with experimental data. 
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Gelatin gels have received great attention 1 * because 
of the interesting fundamental aspects of their rheologi- 
cal behavior and their numerous industrial applications. 
When a semi-diluted gelatin solution is cooled below 
room temperature, chains start to form triple helices and 
progressively a connected network is built. By adding 
crosslinker molecules to the gelatin solution above the 
temperature of the coil-hclix transition, the thermore- 
versible gelation of helices can be avoided and a chem- 
ical gel network is formed, as chain aminoacids react 
with crosslinkers. This possibility of chemical gelation 
in gelatin solutions is extremely relevant for the numer- 
ous applications in pharmaceutical, photographic and 
biomedical industries. The physical gel is characterized 
by an extreme biodiversity due to chemical composition 
of the native collagen, molecular weight distribution, so- 
lution properties such as concentration or pH, which may 
affect the temperature of helix formation^. In addition, 
the option of chemical gelation or a combination of chem- 
ical and physical gelation represent an extremely inter- 
esting way of enhancing the performance of these ma- 
terials, in terms of structural and rheological features. 
As a consequence, an important issue for the develop- 
ment of technological applications is the understanding 
and the governance of the interplay of the different con- 
trol parameters (gelatin concentration, crosslinker con- 
centration, pH or temperature). 

Recently, extended studies have been performed on 
gelatin in solution with bisvinylsulphonemethyl (BVSM) 
reactant 3 , able to establish bifunctional covalent bonds 
with the lysine, the hydroxylysine and possibly with 
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other amine groups of gelatin chains. The influence on 
bond formation of various parameters, i.e. the concen- 
tration of gelatin and reagent as well as the solution pH, 
was investigated. Microcalorimetry measurements were 
able to monitor the development of the chemical reac- 
tion in time by detecting the exothermic enthalpy change 
during the formation of C — N bonds. The kinetics of 
cross-link formation was found to follow a double expo- 
nential decay with two characteristic times. However, 
when counting the number of cross-links binding BVSM 
and gelatin, the method could not discriminate between 
bonds established by free reactants with a chain, bonds 
leading to a bridge between two gelatin chains, or bonds 
leading to a loop within a chain. This lack of informa- 
tion on the kinetics clearly hindered the characterization 
of the structure and therefore the understanding of the 
mechanical properties of the gel, requiring an alternative 
approach to reach a deeper comprehension. A joint in- 
vestigation, pursued by complementing the experimental 
study with suitably designed numerical simulations, al- 
lowed to shed some light on these unanswered questions. 

By means of Monte Carlo simulations on the cubic lat- 
tice of a simple model we analyzed the kinetics of bond 
formation in chemical cross- linking of a gelatin solution 4 . 
We considered a solution of polymer chains at different 
concentrations where reactant monomers can diffuse and 
form bonds with the active sites along the chains, produc- 
ing the cross-linking. It was possible to follow the kinetics 
of the gel formation by varying the gelatin concentration, 
the cross-linker concentration and its bonding probabil- 
ity (i.e. reactivity). The numerical data were able to 
reproduce extremely well the experimental findings. The 
combined analysis showed that the two time scales de- 
tected in the experiments correspond to the average times 
for forming, respectively, single bonds reactant-chain and 
bridges chain-chain via cross- linkers. We related these 
two times to the characteristic times of diffusion of free 
reactants and reactants which have already formed one 
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bond with a chain. Their ratio controls the kinetics of 
bond formation: varying the concentration, the cross- 
linker reactivity and the pH strongly affects this ratio 
and therefore the kinetics of the gelation process. We 
could also show that the reaction rate for a reactant to 
form a bridge between two active sites allows to finely 
tune the kinetics of gelation by affecting the ratio of the 
two characteristic times. This can be achieved in exper- 
iments, for instance, by changing the reactant agent in 
the gelatin solution. 

Another crucial question put forward by the experi- 
mental study^ is the role of the pH. In particular, a qual- 
itative change of the gelation kinetics was observed by 
changing the pH, which could not be explained, also due 
to the limited access to structural information in the ex- 
periments: The ratio of the two characteristic times ex- 
perimentally appears to increase with the solution pH. 
In fact, a higher pH activates more amine groups able 
to react with BSVM along the gelatin chains. This can 
be mimicked in the simulation by varying the number of 
active sites along the chains. Preliminary data obtained 
in simulations^ showed some inconsistencies with the ex- 
perimental observations for gelatin solutions with BSVM. 
Indeed, the possibility of detecting different regimes and 
varying meaningfully the number of active sites along the 
chains strongly relies on the use of long enough chains. In 
this paper we present therefore a new detailed numerical 
study of the kinetics of the bond formation as a func- 
tion of the pH and of the concentration of gelatin and 
reactant. Here longer chains and larger system sizes are 
considered, with the specific aim of modeling solutions 
with different pH and exploring more extended regions 
of the phase diagram, with respect to the previous study. 
We are able to show that the gelation line exhibits a 
non-monotonic behavior, that is, for high reactant con- 
centrations the system goes back into the sol phase for all 
pH values. This result is of course extremely relevant for 
applications and is indeed in agreement with experimen- 
tal findings, where gelation is not observed even at rela- 
tively high cross-linker concentrations^. We show that 
this feature is due to an enhanced probability of saturat- 
ing all the active sites by a single bond with a different 
crosslinker. In specific conditions of high crosslinker con- 
centrations, this can actually hinder gelation, resulting in 
a reentrant phase diagram. By systematically studying 
the dependence of the ratio of the two times on the gelatin 
concentration for different pH, we are able to reproduce 
and rationalize the experimental data. These indeed cor- 
respond to a range of cross-linker reactivity lower than 
the one investigated in the previous study. The present 
results lead to a coherent scenario for the gelation ki- 
netics, which exhaustively complements the experimental 
study. 

The paper is organized as follows: In section U the 
numerical model and the numerical simulations are de- 
scribed in detail. The results on the gelation phase dia- 
gram are discussed in section|TTl where also the number of 
loops and single bonds formed at the end of the reaction 



is analyzed, together with the size of the macromolecule. 
This last feature is not accessible experimentally and is 
clearly extremely relevant for the mechanical properties 
of the gel. In section HTT1 we study the kinetics via the ra- 
tio of the two characteristic times as a function of gelatin 
and reactant concentrations, as well as of the solution pH 
and reactant reactivity. Concluding remarks are given in 
section [IVJ 



I. MODEL AND NUMERICAL STUDY 

We perform Monte Carlo simulations on a cubic lat- 
tice of a system made of bi-functional monomers. Each 
monomer represents a unit on the chain or else a reac- 
tant. The chains are formed by a sequence of n = 20 
linked monomers. One monomer of the chain models 
a Kuhn segment^, and therefore represents more units. 
The length of a Kuhn segment in a gelatin chain has been 
measured 8 to be of the order of 40 A, corresponding to 
about 10 amino-acids. As compared to the experiments, 
our chains correspond to shorter gelatin chains, contain- 
ing only about 200 amino-acids. Each monomer occupies 
simultaneously the eight sites of the lattice elementary 
cell and, to take into account excluded volume interac- 
tions, two occupied cells cannot have any site in common. 
A fixed number of monomers along the chain are active 
sites which may bind to the reactant in order to form 
complex clusters of chains leading to the formation of a 
gel. The active sites are tetra-functional: Two bonds are 
formed with the neighbors along the chain and two are 
not saturated at the beginning of the simulation. The 
number of active sites per chain, n asi corresponds to a 
fixed pH of the solution. Although the number of amine 
groups in a gelatin chain actually linked to reactant can- 
not be measured experimentally, it is estimated that at 
most a fraction of 20% can react. Therefore we perform 
simulations for n as varying from 5 to 11, which corre- 
sponds to a fraction of about 5% to 10% active amino- 
acids in our chain. 

Chains are randomly distributed on the lattice and 
diffuse via random local movements using the bond- 
fluctuation dynamics^. After equilibration, we add the 
reactant to the system and let the solution diffuse to- 
wards the stationary state. Due to the diffusion of cross- 
linkers and chains, when a free reactant finds a nearest 
neighbor unsaturated active site, the first bond forms 
along lattice directions. The second bond is instead 
formed with probability < 1 , since its formation might 
require to overcome a free energy barrier^i^, depending 
on the nature of the solution, the active sites and the re- 
actant. In particular, in these systems the rigidity of the 
C — N link may significantly limit the angle between two 
bonds or the effective chemical reactivity of the cross- 
linker may undergo meaningful variation once the first 
bond is formed. The process goes on until all the possi- 
ble bonds are formed. 

We perform numerical simulations of the model for 



lattice size L = 200, where the unit length is the lat- 
tice spacing a = 1, with periodic boundary conditions. 
The chain concentration C and the cross-linker concen- 
tration C r are defined as the ratio between the num- 
ber of present monomers and the maximum number of 
monomers N max = L 3 /8 in the system. Using the per- 
colation approach we identify the gel phase as the state 
in which there is a percolating cluster—. For a fixed set 
of parameters we generate a number of configurations of 
the system and monitor the reaction. In order to locate 
the gelation transition we analyze the percolation proba- 
bility II, defined as the fraction of configurations leading 
to a percolating cluster, and we identify the transition 
with the line II = 0.5±2. 

We investigate the behavior of the number of bonds 
formed during the reaction process and we classify the 
bonds in three categories: 

1. Bonds between a free reactant and an active site (we 
will refer to this type of bond as single-bonds); 

2. Bonds between a linked reactant and an active site of 
another chain (bridges)] 

3. Bonds between a reactant and two active sites of the 
same chain (which in the following we will call loops). 
We analyze the kinetics of bond formation varying 1 
chain concentration C, the cross-linker concentration ( 
the probability pb of second bond formation and the nu 
ber of active sites n as . The time is measured in Moi 
Carlo unit time, i.e. the attempt to move all monom 
once. The number of bonds is normalized by the ma 
mum number of bonds that can be formed depending 
the limiting agent of the simulation: The normalizati 
factor is the minimum value between twice the number 
reactant monomers and twice the number of active si 

in the system. 



II. PHASE DIAGRAM AND LOOPS 

We first determine a qualitative phase diagram 
varying the chain and cross-linker concentrations, C a 
C r respectively, for n as — 5,8,11, which correspond 
different values of the pH. The results are shown in Fig 

Following the standard site-bond percolation picture- 
one expects that the gel phase is obtained provided that 
enough crosslinkers are added. However, it has been ob- 
served in experiments that this is not always the case, 
and that this feature is affected by the value of the pH. 
FigfT] shows that for all values of n as gelation can be 
obtained in a limited range of C r . Indeed, the figure 
clearly indicates that, at a given gelatin concentration 
C, increasing C r beyond a certain value corresponds to 
a strong decrease of percolation probability and hinders 
the presence of a gel at the end of the reaction. Of course 
this reentrant phase diagram is the product of the bal- 
ance among gelatin concentration, crosslinker concentra- 
tion and number of active sites per chain. Our reen- 
trant phase diagram is in agreement with experimental 
findings for polymer aqueous solutions^, where gelation 
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FIG. 1: The phase diagram, obtained plotting the percolation 
probability II as a function of chain and reactant concentra- 
tion, C and C r respectively, for C ranging from C = 0.01 
to C = 0.14 and CV ranging from CV = 0.001 to CV = 0.14. 
The spanning probability has been averaged over 32 indepen- 
dent initial configurations of a sample of size L — 200 for 
n as = 5, 8, 11 and pb = 1. 




FIG. 2: Normalized number of single bonds reactant-chain 
n s (a), loops m (b) and bridges nbr (c) as function of gelatin 
concentration for n as = 5, pt = 1 and different values of CV. 



can be inhibited by different mechanisms, namely poly- 
electrolyte effects due to the complexation process. 

If the number of active sites per chain increases, n as — 
5, 8, 11, i.e. the pH increases, the gel phase moves towards 
regions with higher concentration of cross- linkers. This 
suggests that in all cases, for sufficiently high C r , the 
system goes back into the sol phase. 

To better understand these results, we measure the 
number of unsaturated reactants, the number of loops 



and bridges at the end of the reaction. Such informa- 
tion, not accessible experimentally, is extremely relevant 
to the mechanical properties of the gel. In Figf2]we plot 
the normalized number of single bonds n s , corresponding 
to the number of unsaturated reactants, the normalized 
number of loops ni and the normalized number of bridges 
ribr (defined at the end of section [TJ) , as a function of the 
chain concentration C for different values of C r . 

We can detect three different regimes: i) At low C r , 
when the system is in the sol phase, the majority of bonds 
are loops at low C and bridges at high C; ii) in the inter- 
mediate regime the number of loops decreases in favor of 
bridges leading to gelation; iii) at high C r the majority 
of bonds are single bonds crosslinker-chain, which brings 
the system back into the sol phase. 

As a small remark, we notice that at very low C r , cor- 
responding to the sol phase, ni decreases with C over 
the whole observed range following a power-law decay 
ni ~ C~ ' 8 , in agreement with previous results 4 . For 
increasing C r , the same power law decay is recovered at 
higher C. The analysis of ni in systems with different pH 
indicates that a higher fraction of loops is measured at a 
higher pH, as expected. 




FIG. 3: Fraction of monomers in the largest cluster for n as — 
5 and different reactant concentrations, ranging from C r = 
0.001 to C r = 0.1. The symbols indicate the same values of 
C r as in Figg] 

We also measure the fraction of monomers belonging 
to the largest cluster, p max = s max /N, where N = (C + 
C r )L 3 /8 is the total number of monomers in the system. 
Here s m ax is evaluated by counting monomers, either in a 
chain or reactants. In Fig\3\p max is plotted as a function 
of C. Except for the very small C r , corresponding to 
the sol phase, the fraction of monomers belonging to the 
largest cluster increases toward unity for increasing C, 
as the system moves into the gel phase. For higher C, 
Pmax exhibits the tendency to decrease, as the system 
goes back into the sol phase. 

Discussion: Depending on specific conditions, the 
complex balance between C, C' r and n as , results in en- 
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FIG. 4: (Number of single-bonds (a), number of bonds form- 
ing bridges or loops (b) and total number of bonds (c), nor- 
malized by the total number of possible bonds, as a function 
of time for C = 0.1, CV = 0.025 with n aB = 5 (circles), 8 (tri- 
angles), 11 (squares) and pi, = 1. At the end of the reaction 
all systems are in the gel phase. 



hancing the probability of forming either loops or sin- 
gle bonds reactant-chain. For all pH values and low C r 
gelation is not observed as reactants mostly form loops 
within the same chain. Conversely, gelation is hindered 
at high C r since reactants saturate all the active sites 
of the chains, not leading to the formation of a macro- 
molecule. Indeed this effect is observed at higher C r for 
higher n asi since, at higher pH, more reactants arc needed 
to saturate all active sites. Therefore, beyond a certain 
value of C r , adding reactants to the solution is useless, 
as a higher concentration of chains is needed for gelation. 
In the following section we investigate how the complex 
scenario above described is connected to the gelation ki- 
netics. 



III. KINETICS OF BOND FORMATION 

In experiments the pH of the solution is found to 
strongly influence not only the final state of the system, 
i.e. the phase diagram, but also the kinetics of bond 
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FIG. 5: Number of single-bonds (a), number of bonds forming 
bridges or loops (b) and total number of bonds (c) , normalized 
by the total number of possible bonds, as a function of time 
for C = 0.01, C r = 0.025 with n as — 5 (circles), 8 (triangles), 
11 (squares) and pb — 1. At the end of the reaction all systems 
are in the sol phase. 
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FIG. 6: Ratio between the two characteristic times for n as = 
5,8,11, p b = 1 and for CV = 0.001 (a) CV = 0.05 (b) and 
C r — 0.05 (c). The arrows indicate the value of gelatin con- 
centration where gelation occurs for the different numbers of 
active sites. 



formation. 

In order to investigate this phenomenon, during the 
reaction we monitor the total number of bonds formed 
nb(t), the number of single-bonds n s (f) (bonds of type 1) 
and the number of bonds forming bridges or loops n<i (f ) 
(i.e. bonds of type 2 and 3), with obviously n b (t) = 
n s (t) + ri2 (t) . The number of bonds is normalized by the 
total number of possible bonds. 

In Figf?]n s (t), 722(f) and n&(f) are plotted as a function 
of time in the case pb = 1 for n as = 5,8, 11 for C = 0.1 
and CV = 0.025, corresponding to the gel phase at the 
end of reaction. As the reaction begins, single-bonds form 
rapidly, then second bonds start to form and the degree of 
connectivity between chains increases: single-bonds de- 
crease in time as the number of bridges increases. Our 
data indicate that as n as increases, second bonds form 
faster, whereas the characteristic time of single-bond for- 
mation is less affected by pH. As a consequence, we ex- 
pect that the ratio of the two times decreases with pH. 
A significantly different behavior is observed in the sol 
phase (FigfSJ): in this case the characteristic time for the 
first and the second bond formation are the same for all 



pH values. Indeed, in this regime the whole kinetics of 
bond formation is independent of the solution pH (see 
FigJSJ:). Hence, for pb = 1, the time ratio is expected to 
be almost constant by varying the pH. 

We now focus on the ratio between the characteristic 
times n of single bond formation and T2 of second bond 
formation and we study its dependence as a function of 
the number of active sites per chain n as . In FigJS] we 
distinguish different regimes, according to the concentra- 
tion of chains and cross- linkers. For low C r (FigJBfa)) all 
systems are in the sol phase due to the very high number 
of loops (Fig(2J) and therefore the time ratio is close to 
one. As the number of active sites per chain increases, 
Ti remains quite constant whereas T2 decreases, so that 
T%jT\ slightly decreases as n as increases. For intermedi- 
ate cross-linker concentrations (Fig^b)) and low chain 
concentrations (sol phase), neither t\ nor T2 strongly 
vary increasing n as , consequently the ratio Tijr\ is al- 
most independent of n as and weakly decreases increasing 
n as . For higher chain concentrations, beyond the gela- 
tion threshold (indicated by an arrow in the figure), the 
time ratio starts to increase, due to the slower decrease 
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FIG. 7: Ratio between the two characteristic times for 
n as = 5, C r = 0.025 and p b — 0.001 as a function of chain 
concentration C. 



of t 2 with respect to n, caused by the presence of the 
macromolecule. The smaller the number of active sites 
per chain, the later the increase of r 2 /t± is observed with 
the chain concentration, in agreement with the displace- 
ment of the gel line in the phase diagram of FigfT] More- 
over, for smaller number of active sites the height of the 
peak decreases, suggesting that r 2 decreases faster than 
T\. For high C r (FigJS^c)) the position of the maximum 
time ratio appears to move toward higher chain concen- 
trations, due to the shifting of the gelation transition. 

We conclude our study by investigating the interplay 
between the cross-linker reactivity and the pH. A first im- 
portant hint of the possible combined effects is obtained 
by comparing the data of Fig|6jb), discussed above and 
corresponding to Pb = 1, to the behavior of t^/tx as a 
function of C measured for the same C r = 0.025 but 
for a lower reactivity reactant, pb — 0.001 (Fig[7]). In 
this case the time ratio is about one order of magnitude 
larger than for high reactivity crosslinkers (Fig[S]). More- 
over, the dependence of Tijrx on the number of active 
sites seems to be less strong and qualitatively different. 
In fact, in contrast to what observed in FigJBJb), here in 
the gel phase t 2 /ti increases with increasing n as , since 
Ti decreases whereas r 2 is very large for all n as due to 
the very small pb- 

In general, if the probability pb of forming the second 
bond varies, the mean time of second bond formation r 2 
changes, and so does the velocity of the reaction. In Fig[5] 
we plot the ratio t 2 /ti as a function of pb for different 
values of n as . It is worth to notice that the average n 
is constant along each line, since pb only affects t 2 . If 
the bond probability is very low (p& < 0.005) an increase 
of n as will correspond to an increase of r 2 /t\ , in agree- 
ment with Fig. 7. In this case the ratio is controlled by n, 
which decreases as the number of active sites increases. 
Conversely, if pb is sufficiently high (pb > 0.01) we ob- 
serve the opposite behavior (see Fig. 6(b)). In this case, 
the ratio is controlled by r 2 : As the number of active 
sites per chain increases, the mean time to form a second 
bond r 2 decreases more rapidly, affecting the ratio r 2 /n . 
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FIG. 8: Ratio between the two characteristic times for n aa = 
5,8, 11, C r = 0.025 and C = 0.1 as a function of bond prob- 
ability pt,. For all n aa the system is in the gel phase. Inset: 
Ratio between the two characteristic times as a function of 
solution pH in of gelatin chains in solution with BVSM. 



Discussion: Our data indicate that when the sys- 
tem at the end of reaction is in the sol phase, the ratio 
between the characteristic times does not substantially 
vary with the pH. Conversely, in the gel phase the kinet- 
ics of bond formation is strongly influenced by the n as , 
i.e. by the pH. Since the time ratio t 2 /ti is actually mea- 
surable from the experiments^, our observation suggests 
a new alternative way to discriminate the sol from the 
gel phase resulting from crosslinking the gelatin solution 
in specific conditions of C and C r , purely based on the 
reaction kinetics. 

It is worth to remark that the presence of loops, in 
general, corresponds to smaller values of t 2 , whereas the 
formation of bridges leads to larger values, as follows from 
the microscopic interpretation of t\ and t<£. At high C 
the effect of the pH in very dense systems becomes less 
important. Indeed, at high C both t\ and r 2 decrease 
as n as increases, so that the ratio t 2 /ti tends towards a 
value which is independent of the number of active sites 
per chain. This analysis allows to identify the range of 
parameters C and C r for which the reaction kinetics is 
more sensitive to pH changes. This is certainly an impor- 
tant information for the experimental characterization of 
these systems. Moreover, it also suggests possible new 
applications in terms of pH-responsive systems. 

Furthermore, for low reactivity cross-linkers the kinet- 
ics of the reaction is mainly controlled by the charac- 
teristic time for forming the first bond, and therefore 
the time ratio increases with the pH. For high reactiv- 
ity cross-linkers, conversely, the combined effect of the 
pH and of the reactivity makes r 2 decrease faster than 
Ti and therefore a decrease of the ratio will be observed 
when n as increases (i.e. for higher pH). By measuring 
the dependence of the time ratio on the pH, it is then 
possible to infer the level of reactivity of the reactant. 
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For instance, experimental data for gelatin solutions with 
BVSM reactant^ show that, at fixed C and C r , the time 
ratio increases with the pH (Inset Fig|8|). According to 
our analysis this would suggest that BSVM is a low re- 
activity cross-linker, namely that the number of amine 
groups and reactant configurations leading to a chemi- 
cal bond is small. Interestingly enough, this finding does 
complete an independent previous analysis^, which sug- 
gested that the peculiar two-time kinetics observed in 
these systems should correspond to a low value of the 
crosslinker reactivity pi, . 

IV. CONCLUSIONS 

In conclusion, we have performed a comprehensive 
study of the pH dependence of the phase diagram and the 
reaction kinetics in crosslinking gelatin solutions. This 
study is based on numerical simulations of a lattice model 
specifically developed for these systems. The results al- 
low to rationalize the experimental findings and give new 
relevant insights. The first important result of our study 
is that the increase of the cross-linker concentration does 
not necessarily imply the presence of a gel at the end of 
reaction: indeed it may also hinder the gelation transi- 
tion, inducing the formation of a high number of single 
bonds or loops. This effect moves towards higher con- 
centrations if the pH of the solution increases. As a con- 



sequence, we observe a reentrant phase diagram. 

Our study points out that also the kinetics of bond for- 
mation is strongly influenced by the pH of the solution, 
with different regimes according to the concentrations of 
chains and cross- linkers. We find a remarkable qualita- 
tive difference in the pH dependence of the reaction kinet- 
ics, which could be used as a new purely kinetic criterion 
to discriminate between the sol and the gel final state of 
the crosslinking process. Thanks to a detailed analysis 
of reaction kinetics in terms of single bonds, loops within 
chains and bridges between different gelatin chains, we 
are able to give a meaningful explanation of the different 
regimes observed. Most interestingly, our results show 
that there exist ranges of the control parameters C and 
C r for which the behavior of the system may be extremely 
sensitive to pH changes. These findings are fundamen- 
tally important for the experimental studies and suggest 
new possible applications of these materials. 

Finally, we have been able to clarify how the influence 
of the pH on the kinetics of bond formation also depends 
on the reactivity of cross- linker considered, i.e., on the 
bond probability p^. For low reactivity cross- linkers the 
kinetics of the reaction is mainly controlled by the single- 
bond formation time, which decreases with pH, leading 
to an increase in the time ratio. Conversely, for high 
reactivity cross-linkers and for high chain concentration, 
the increase of pH induces a definite decreases in the time 
ratio. 
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